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Abstract. We further develop the Tayler-Spruit dynamo theory, based on the most efficient instability for gener- 
ating magnetic fields in radiative layers of differentially rotating stars. We avoid the simplifying assumptions that 
either the fi- or the T-gradient dominates, but we treat the general case and we also account for the nonadiabatic 
effects, which favour the growth of the magnetic field. The general equation leads to the same analytical solutions 
in the limiting cases considered by Spruit ( 2002J. Numerical models of a 15 Mq star with a magnetic field are 
performed. The differences between the asymptotic solutions and the general solution demonstrate the need to 
use the general solution. Stars with a magnetic field rotate almost as a solid body. Several of their properties (size 
of the core, MS lifetimes, tracks, abundances) are closer to those of models without rotation than with rotation 
only. In particular, the observed N/C or N/H excesses in OB stars are better explained by our previous models 
with rotation only than by the present models with magnetic fields that predict no nitrogen excesses. 
We show that there is a complex feedback loop between the magnetic instability and the thermal instability driving 
meridional circulation. Equilibrium of the loop, with a small amount of differential rotation, can be reached when 
the velocity t/magn of the growth of the magnetic instability is of the same order as the velocity C/circ of the 
meridional circulation. This opens the possibility for further magnetic models, but at this stage we do not know 
the relative importance of the magnetic fields due to the Tayler instability in stellar interiors. 
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K*" 1. Introduction is the main restoring force, i.e., when A'^, < N^. In this 

case, the radiative losses from the magnetic instability are 

A major concern is to know whether magnetic fields are ^j^^ accounted for; these losses reduce the restoring force 

^ . important for stellar evolution or not (cf. Roxburgh 2003) . buoyancy which acts against the magnetic instability. 
■ ■ ' Some initial numerical tests have suggested that the mag- 

netic fields created by the Tayler-Spruit dynamo (Spruit The problem with the above description is that, in the 

[TM9. 20021 in rotating stars may be large enough to trans- interiors of models no longer on the zero-age sequence, 

port angular momentum in stellar interiors as well as ^^^re is at some moderate distance of the core a signif- 

chemical elements in upper MS stars (Maeder & Meynet i^^nt zone with < Nt (covering about 25 % of the 

mm . In this work, we perform further theoretical devel- stellar radius), which is in a situation where the current 

opments of the dynamo process, as well as some numerical simplifications made in Case 1 are not appropriate. In this 

models, to study the consequences of stellar evolution with ^o^e (we called it zone IP), the magnetic instability is 

a magnetic field adiabatic, while the hypothesis of large non-adiabaticity 

T • .1, r i-- ■ i- un-4-- j-i is made (see Sect. 2 below). It also happens in all stars 

In the nice theory of magnetic instabilities recently ^ ' i r i 

, , 1 1, c «Tw™ 1 • J- J- T v that there are zones where Nn and Nt are of the same 

developed by Spruit (|2002(l . only some important limiting ^ i i i 

,1 -J J T 1 j-T, 11 J order and where it is better to account for both gradients, 

cases have been considered. In particular, the so-called , , , „ . , ' 

^ „ , ii-iT/-i even it m general such zones are rather small. We consider 

Cases and 1 have been studied, in Case 0, the /^-gradient . ^ , , . , , , , 

J . , 1 J- i • i 11 • i • it preferable to consistently treat the general case and to 

dominates over the thermal gradient m stellar interiors, , r i i 

. , r , , 1 TD J- r- r account properly everywhere for the T- and u-gradients, 

I.e., m terms of the partial Brunt-Vaisala frequency we u r i ii rr 

u AT ^ T\T 1 1- 1, -tu 1 J- 4. as well as for the non-adiabatic effects, 
have Nfj^ > Nt- Case 1 applies when the thermal gradient 

Also, the two approximations mentioned above have 
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lar parameters, such as the angular velocity Q, the differ- 



ential rotation q = 



dinn 

J In r 



or the thermal and chemical 



gradient. The use of asymptotic solutions such as Spruit's 
H2(JU2|I introduces artificial jumps in the behaviour of the 
diffusion coefficients. 

Our purpose here is to express the theory in a unique 
and consistent formulation for all possible cases of Nfj_ and 
A^T and of departures from adiabaticity. In Sect. 2, we dis- 
cuss the non-adiabaticity as a function of radiative and 
magnetic diffusivities K and 77 respectively. In Sect. 3, 
we establish the general equations for the magnetic field 
and we check for consistency with the expressions derived 
by Spruit (j2002|l . In Sect. 4 we examine the general ex- 
pressions for the transport of the angular momentum and 
of the chemical elements. Sect. 5 shows some initial nu- 
merical applications as an example. In Sect. 6 we exam- 
ine the reciprocal feedback between meridional circulation 
and Tayler instability and the fact that these two processes 
lead to very different amounts of differential rotation. The 
conclusions are given in Section 7. 

2. The radiative heat losses of the magnetic 
instability 

2.1. Magnetic difFusivity and radiative losses 



As shown by Spruit H2()()2|l . the magnetic diffusivity 77 in 
Case (A^^ > A'^t) is given by 



riQ = r n 



while in Case 1 {Nf^ < iVr), it is given by 

The Alfvcn frequency wa is given by 
B 



(1) 



(2) 



(3) 



The Brunt -Vaisala frequency, when the /x-gradient dom- 
inates over the T-gradient, tends towards the following 
partial frequency 



^ Hp ^ ' 



(4) 



and when the T-gradient dominates, it tends towards 



Nr 



Hp 



(Vad - V) 



(5) 



The thermodynamic coefficients S and tp are defined as 

follows: 5 = and ^ = {§^)^^^- Hp is 

the pressure scale height. We need a general expression 
for the magnetic diffusivity encompassing the two limit- 
ing cases with 770 and 771. This requires further study of 
the nonadiabatic heat losses. For a fluid element displaced 



over a length I, the thermal fluctuations diffuse with a 

.2 

timescale ^therm ~ 7^ , where K is the thermal diffusivity 



K = 



Let us call t the typical timescale over which 
the magnetic instability develops. In a time t, the fluctua- 
tions are reduced by a factor /, which can be written (cf. 

Spruit mna), 



/ - 



^thc 



1 



(6) 



If / ^ 1 as for <thcrm ^ t, whcrc t is the dynamical 
timescale of the instability, one has the adiabatic case. On 
the other hand, if / — > 00 as for ttherm *C t, we have a 
highly nonadiabatic case. Spruit H2()()2|l also introduces in 
Case 1 (when A^t dominates) an "effective thermal buoy- 
ancy frequency" A^c = ^ and assumes ithorm ^ t. This 

implies that / ~ 7—^ — . If it happens in some stellar layer 

'■therm 

that t <C ithcrm: situation is adiabatic and one should 
have / = 1, but the adopted simplification leads to / = 0. 
The approximation is thus invalid in this case. 

The factor / introduced by Spruit is related to the 
usual Peclet number Pc = ~ which is the ra- 

tio of the thermal to the dynamical timescales for a fluid 
element of velocity v displaced over a length I. In astro- 
physics, rather than the Peclet number, one often consid- 
ers (e.g. in the theory of convection) a factor F, which is 
the ratio between the energy delivered by a fluid element 
to the energy lost during the displacement of the fluid el- 
ement. For a typical geometry, i.e., by assuming that the 
fluid elements are spherical, one would have F = ^ Pg. 
The fluid elements are likely not spherical in the relevant 
geometry for Tayler instability. We should rather consider 
thin slabs of thickness I, (we are indebted to Henk Spruit 
for this remark). If so, the relation between F and the 
Peclet number is F = i Po- Thus, the following relation is 
valid for the factor / introduced by Spruit and the usual 
number ratios Pn or F 



F.i- 
2 (/ 



1 



or / = 



1 



2 F 



2 F 



(7) 



This comes directly from the abovementioned definition of 
the Pcclct number as well as from Eq. ® . In the adiabatic 
case, / — > 1 and F —^ 00; if non-adiabatic effects are 
important, / is large and F is very small. 

2.2. Oscillations frequencies 

The general expression of the Brunt- Vaisala frequency for 
a displaced fluid element is 



9l 
Hp 



(V - V 



(8) 



where V' is the internal gradient of the the fluid ele- 
ment, i.e., expressing how T changes inside the fluid ele- 
ment during its motion, while V is the T-gradient of the 



Andre Maeder, Georges Meynet: Stellar evolution with rotation and magnetic fields: 



3 



surrounding medium. One also has the following relations 
between the various T-gradients (of. Maeder [l995|l 



V- V = 



r + 1 



(Vad - V) 



(9) 



and 



V = 



V 



rad 



1+r 



'ad 



1 



2 

i+r 



with a numerical factor of 6 instead of 2, if we would have 
considered spherical fluid elements. Thus, in a radiative 
zone where there is some motion of matter, which may 
be non-adiabatic, the thermal gradient depends on the 
eSiciency F of the thermal transport and the gradient may 
lie anywhere between the adiabatic and radiative gradients 
depending on this efficiency. If we write T in terms of /, 
we have 



r 



1 



1 



2/-1 



3. General equation for the magnetic difTusivity 



As emphasized by Spruit f2Qil2j) , the energy of the Tayler 
instability (Tayler J.973 ) must be large enough to overcome 
the restoring force of buoyancy and this implies that the 
Alfven frequency must be larger than some limit. Thus, 
the vertical extent of the magnetic instability must be 
limited (cf. Spruit ESEl Eq. 6), 



(10) k < 



r LP A 
N 



(15) 



However, if this radial scalelength is too small, the per- 
turbations will be quickly damped by the diffusion of the 
magnetic field characterized by i] and this implies > ^ , 
where ctb is the characteristic frequency of the magnetic 
field. In the absence of Coriolis force, this frequency is 
(Jb = ^a- However, as shown by Spruit H2UU2I see also 
Pitts & Tavler [T986|l . in a rotating star (typically when 
> oja) this frequency becomes ctb = due to the 

Coriolis force. Therefore one has in this case 



(11) 



> 



O-B 



(16) 



Thus, we may also write the Brunt-Vaisala frequency tak- 
ing account of radiative losses in terms of /, as follows 



1 



2/- 1 



(12) 



When A^^ = 0, Spruit ij^nO^ writes N'^ which has a 

dependence on 1// for large values of /; however we note 
that the numerical factors may introduce some differences. 

We can go a step further in expressing T in term of 
the ratio "q/ K of the magnetic to the thermal diffusivity. 
As seen above, one has F ~ i where I is the length 
scale of the instability. Instabilities over smaller length- 
scales are removed by magnetic diffusivity. We consider 
here the marginal case P = rjt (cf. Spruit I7U(J2|I and thus 
we get 



1 T] 

2 H 



Thus, one obtains the ratio 



(13) 



which leads to 



the following expression for the Brunt -Vaisala frequency 



N"^ = ^!-^ ^ + A^^ 



T]/K + 2 



(14) 



This is the oscillation frequency of a fluid element dis- 
placed by Tayler-Spruit instability. The numerical factor 
of 2 in the above expression is due to the geometry we 
considered; in many cases this numerical factor will intro- 
duce some minor differences with to Spruit H2002|l . The 
resulting differences will be rather small because this fac- 
tor will generally appear with a power smaller than 1, as 
for example in Eq. 135(1 . 



The combination of these two limits together with Eq. H14|l 

leads to the following condition {^)* > ^ WTi^ where 
N"^ is given by Eq. (|14|l . For the marginal situation corre- 
sponding to equality, we have 



A^ 



(17) 



This equation relates the two unknown quantities, the 
magnetic diffusivity rj and the Alfven frequency uia- 

For the particular cases studied by Spruit H2l)U2|l . the 
result was a much simpler equation. Here, we find it prefer- 
able to establish a second relation between the two un- 
known quantities. The amplification time Ta of the Tayler- 
Spruit instability, i.e., "the timescale on which the radial 
field is amplified into an azimuthal field of the same 
order as the already existing azimuthal field" can be writ- 
ten Ta = N/{ujAVLq) (Spruit JO(52|). The equality of the two 
timescales Ta and Tg^, according to the expression given 
above, leads to a second equation ^ = g ^, as found by 
Spruit H2UU2I his Eq. 18). When account is taken of the 
expression of the Brunt- Vaisala fEa. I14|l . this gives 



A2 



n/K 



(18) 



The above two equations (|17II18|I form a coupled system of 
degree 4 (see below) with two unknown quantities, rj and 
ujA- This system could be solved by standard procedures, 
however we first make some simplifications. We eliminate 
the expression of A^^ between these two equations and 
obtain 



72 Kn ) 



(19) 
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This is a new general equation connecting the mag- 
netic diffusivity and the Alfven frequency. It shows that 
the Alfven frequency, and therefore the magnetic field, 
changes as 1/6 the power of the magnetic diffusivity for 
a given rotation. This coefficient applies to the diffusion 
of the azimuthal component of the magnetic field. Due 
to the particular pinch-type nature of the Tayler-Spruit 
instability, it also applies to the turbulent diffusive mix- 
ing of the chemical elements by this instability. To ensure 
consistency, we verify that in Spruit's Case with 



/WAV 

\n Jo 

we get 



(20) 



(21) 



This is in agreement with Eq. (42) by Spruit l|2UU2|l . This 
equation shows that the mixing of chemical elements de- 
creases very strongly if the //-gradient grows and this ef- 
fect limits the chemical mixing of elements by the Tayler- 
Spruit dynamo in the regions just above the convective 
core. 

In Case 1, when thermal losses are accounted for, the 
Alfven frequency is given by 



K 



(22) 



This is given by Spruit H2UU2I his Eq. 19). Below, we will 
verify that our system of equations also leads to this ex- 
pression (see Eg. 1351 below). For now, we check that our 
Eq. (|19|l with this expression of (77)1 leads to 



m 



K 



(23) 



We have transformed our system of two equations of de- 
gree 4 with two unknown quantities 77 and loa into one 
equation of degree 4 with only one unknown quantity x. 
The solution will provide the value of the Alfven frequency 
and the magnetic diffusivity 77 by Eq. H19() . Higher order 
terms depending on F = might be considered, but as 
also shown by Spruit, this ratio is very small. The above 
equation applies to the general case where both and 
A^T are different from zero and where thermal losses may 
reduce the restoring buoyancy force. This full equation is 
to be solved numerically (see Sect. 5 below). 

Higher order terms depending on F^ could be consid- 
ered, but it is correct to limit the development to the first 
order in F. Indeed, we notice that the expression of TVt 
given by Eq. 101 contains the term V, which is given by 
Eq. (|10|l . This term V contains the ratio If we would 
fully develop the above equation taking this additional de- 
pendence into account we would get an equation of order 
10 instead of order 4 as above. We remark that the depen- 
dence in introduced by V in the expression of N"^ is of 
an higher order in than the one expressed by Eq. (|14|) . 
Since ranges between 10~^ and lO^'' in the numerical 
models of Sect. 5, we may ignore this higher order term. 

We can also check that Eq. H25() leads to simple 
solutions fully consistent with the results by Spruit 
l)2()()2|l in the particular cases he has studied: 



- 1. Firstly, we notice that if terms going like the power 
6 or higher of (^) are negligible, the solution is simply 



n 



(26) 



which is the same as Eq. (|20|) above. This is the so-called 
Case 0. Indeed, the condition lja ^ ^ leads to a solution 
which is the same as when Nt ~ 0, which defines Case 0. 



which is the same as Eq. (43) by Spruit H2UU2|I . Thus, 
we have checked the consistency of the general expression 
for ri (Eq. with the results obtained by Spruit H2U02|I 
in the two particular cases he considered. There are no 
new physics here with respect to Spruit's work P002), 
we simply verify the consistency of our general expression 
with the asymptotic cases considered by Spruit. 

Now, we can introduce the general expression of 77 
given by (|19|l in the second equation (Eq. I18|) of our cou- 
pled system and obtain 



n4 



Nf. 



q^K \ n J 



2. The general Eq. (|25|l can also be written 



-a;" I —{N4 + N'p-n'j+2q' \^—J^-n^j =0.(27) 

We notice some similarity of the two terms in parentheses. 
If we consider Spruit's Case with Nt — 0, this equation 
becomes the product of two terms 



2 1 ^2 



Since the ratio (^) always appears with a power of 2, we 
now define a new variable x — {^)'^ and thus we get 



q^K 



N4 



K 



-x^ + 2Nl X 



2n^q^ = . (25) 



— x' + 2q' 



= 



(24) 



The solutions are 



iV2 
A' 



and 



x^ = 



2 q^K 



(28) 



(29) 



The second possibility leads to negative solutions, which 
are thus not physically meaningful. Thus, this shows the 
uniqueness of the first solution. This first solution is that of 
Case 0, given above by Eq. (|20|) . as found by Spruit 
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Eq. (EZIl 



3. Let us consider Case 1 with iV^ = 0. We have from 



2 q^n^ = 



If q is small, one has the solution 



(30) 



(31) 



which is the alternative solution proposed by Spruit f2UU2), 
when the simplification resulting from the assumption of 
a large / value does not apply, (this is the case we have 
called IP). 

We may also examine the general solution when the 
ratio fj/K is small, which is generally the case as seen 
above and as mentioned by Spruit. For that, we start from 
the two equations 1171 118|l and make the simplifications 
appropriate for a small ij/K ratio. With developments like 
those done above, we get instead of Eq. (|25(l the following 
expression 



4. Transport of angular momentum and chemical 
elements 

Let us specify some other useful expressions such as the 
critical lengthscale for the Spruit-Tayler instability, the 
intensity of the magnetic field, the coefficients of transport 
for the angular momentum and for the chemical elements. 
As shown above, the vertical extent of the magnetic insta- 
bility is limited by Ir < fEa. [T5|l . where N is given by 
the general expression (|14|l . Thus the maximum length- 
scale of the magnetic instability is 



r ujA 



ri/K+2 



(36) 



Following Spruit H2002|l . we consider this maximal value 
as the one charateristic, for example, of mixing, since the 
mixing will be essentially determined by the instabilities 
with the longest lengthscale. We immediately verify that 
Ir is in agreement with Spruit's H2l)(J2|l result in Case 0, 
when TVt = 0. In Case 1, when = 0, we may write, if 
we assume in addition that is small, 



2^ 



q^K 



:N4 



■2Nlx- 



2 n^q^ 



0, 



(32) 



r UA 



(37) 



which is similar to our Eq. H25|l for iV^ = 0, except that the 
term in has disappeared. Indeed, we may also derive 
the above equation directly from our general Eq. H25|) by 
noting that if rj/K is small, then from Eq. H18|l one has 



and then by multiplying this equation by — ^ , we get 



With the general expression of rj given by Eq. H19|) . this 
becomes 



25 K2 q 



(33) With Eq. (123 for uja, we get finally in Case 1 



, nY f K 



(38) 



(39) 



q^ K 



:N„, 



K 



x-" = 0. 



(34) 



When we bring this relation into Eq. (|25|l . we are lead 
to the above Eq. IpT^ . Thus, this relation (p?^ is fully 
consistent with our general equation. 

In Case 1 when N^, the above Eq. H32|l leads to 



25 



At 



K 



which is identical to the result by Spruit 1)20021 Eq. 10), 
apart from the geometrical factor 2*. Thus, our general 
expression for the lengthscale of the Tayler-Spruit in- 
stability also contains the limiting cases studied by Spruit. 
We also recall, as shown in Eq. (12) by Spruit, that -q <^ K 
is a condition of validity for the developments made by 
Spruit. 

Turning now to the general expression for the intensity 
of the magnetic field, following Spruit H2()()2|l we have for 
the azimuthal and radial field strengths, 



= (47rp)2 r uja and = B^ (k/r) 



(40) 



2« 



At 



K 



r^NT 



(35) 



This is just the same expression, except for the small ge- 
ometrical factor 2s, as in the equation found by Spruit 
ipr02. i.e., his Eq. 19). 

Thus, for all situations considered by Spruit, we have 
verified here that our general equations H191 125(1 well re- 
produce his results. 



The quantity uja has to be taken as the solution of the 
general equation 125|) and ^i- is given by Eq. (|36l) . These 
general expressions also reproduce correctly the particular 
Cases and 1, since we have verified the consistency of our 
general expressions for /r and uja with those for the two 
particular cases. For Case 0, this gives (cf. Spruit H2002|l . 



B^ = iATTp)irq— 



and — - = 

Bip 



A^ 



(41) 
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For Case 1, one has 
B 



K 



r2 Nt 



and 



23 



K 



Nrr 



(42) 



(43) 



These last two expressions are identical to those by Spruit 
except for the geometrical factors 2 8 and 2* respectively. 

Let us now turn towards the transport of angular mo- 
mentum by the magnetic field. The azimuthal stress by 
volume unity due to the magnetic field is given by 



S — - — B^Bu, — - — ( — 

4 TT 4 TT V r 



p r 



Following the ingenious procedure devised by Spruit 
(|2Ul)2|l . we can express the viscosity v for the vertical 
transport of angular momentum in terms of S. Indeed, 
the viscosity of a fiuid represents its ability to transport 
momentum from one place to another and thus we get, 



S _ f7r2 /WA\3 (VL 



(45) 



This is the general expression of v with to a given by the 
solution of Eq. and with N by Eq. . We immedi- 
ately verify that in Case with 77 = 9 we find the 
same as Spruit H2002|l . 



(46) 



The /i-gradient, through its reduction of the field, partic- 
ularly in the vertical direction, also reduces the transport 
of angular momentum, but much less than for the chem- 
ical elements as given by Eq. H21|l . In Case 1, it is also 
easy to verify, by using the appropriate expressions for N, 
and OJA that one obtains. 



Nt 



K 



(47) 



which is identical to Eq. (32) by Spruit H2(J(J2fl apart from 
the geometrical factor 62 . 

Thus, we have the full set of expressions necessary to 
obtain the Alfven frequency wa, the magnetic diffusiv- 
ity ?7, which is necessary for calculating the transport of 
the chemical elements, and the viscosity v for the ver- 
tical transport of the angular momentum by the mag- 
netic field. These expressions apply to the general case, 
where both the /i- and T-gradients are different from 
zero and where radiative losses may play a role. In all 
cases, the general expressions may also lead consistently 
to the asymptotic expressions by Spruit ()2()[)2|l . We may 
remark that we need not to include here the erosion of the 




Fig. 1. Internal distribution of the angular velocity 0(r) 
as a function of the radius in solar units in the model V3 
(without magnetic field) at various stages of the model 
evolution indicated by the central H-content Xc during 
the MS-phase. The initial mass is 15 Mg and Z = 0.02. 

/i-gradient by the horizontal turbulence as suggested by 
Talon & Zahn H1997|l . since as shown by Maeder & Meynet 
(2003), the horizontal turbulence is either suppressed or 
at least strongly reduced by the magnetic field. 

Finally, let us check consistency. The rate of magnetic 
energy production Wb per unit of time and volume must 
be equal to the rate of the dissipation of rotational en- 
ergy by the magnetic viscosity v as given above. This check 
of consistency is verified for the asymptotic expressions, 
but we have to see whether the general expressions also 
fulfill it. We assume here that the whole energy dissipated 
is converted to magnetic energy, which is a reasonable ap- 
proximation. The sink of energy due to chemical mixing is 
negligible here, since as shown below the transport of ele- 
ments from the core to the surface is absent from models 
with the magnetic field as calculated here. 



which gives with Eq. H45|l . 

With Eq. (|18|l defining the field amplitude 
the dissipation rate finally becomes 



(48) 



(49) 



O q 
N ' 



(50) 



We turn to the rate Wb of magnetic energy creation. The 
magnetic energy per unit volume is , it is produced in a 
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Fig. 2. Internal distribution of the angular velocity fl{r) 
as a function of the radius in solar units in the model Ml 
(with magnetic field calculated according to Spruit 2002 ) 
at various stages of the model evolution indicated by the 
central H-content Xc during the MS-phase. The initial 
mass is 15 M,^ and Z = 0.02. 



characteristic time given by 
has 



Wb 



r2 , ,2 

8tt n 



1 



-pr 



2^ 

n 



Lol/n)-'^. Thus, one 



(51) 



where we have used the Eq. (|40() for the field, because 
is the main field component. If we now use the above 
expression of we get the same expression as for W,^ 
fEg.lSnjl. thus one has 



= Wb. 



(52) 



This shows the consistency of the field expression for B^, 
of the transport coefficient v together with the energy con- 
servation in the process of field creation. 

5. Numerical applications 

5.1. Existing models, present models and their aims 

Many sets of models with rotation exist, often without 
magnetic fields. A general review of the models with- 
out magnetic field has been made by Maeder & Meynet 
H2000|l and more recent references are given in the re- 
cent lAU Symposium 215 on Stellar Rotation (Maeder & 
Eenens l^004p . Most models cover the H- and He-burning 
phases, while some models go up to the pre-supernova 
stage (Heger et aL lMTUl Hirschi et al'^UM). The models 
often differ significantly in their physical assumptions and 



Fig. 3. Internal distribution of the angular velocity 0(r) 
as a function of the radius in the model M5 (with magnetic 
field calculated according to expressions of the present pa- 
per) at various stages of the model evolution indicated by 
the central H-content Xc during the MS phase. The dot- 
ted line shows the values at the end of the H-burning 
phase. 



the way they are modelized. A particularly critical point 
is the treatment of the meridional circulation, which in 
many sets of recent models is treated as a diffusion process, 
which is incorrect since a diffusion process goes along the 
gradient of the considered quantity (for example fl) , while 
the advection due to meridional circulation may transport 
the angular momentum from the location where there is 
less (central regions) to the location where there is more 
(superficial regions), and it can also do the opposite de- 
pending on the circulation patterns. 

Several stellar models with rotation and magnetic 
fields have been published in the recent literature. The 
interaction between meridional circulation and dynamo 
action have been studied by Charbonncau & MacGregor 
IpMl see also Sect. 6 below). Hcger et al. ipiUHjl have 
implemented Spruit's H2()()2|l asymptotic developments in 
models of stellar evolution. These authors follow from the 
ZAMS to the pre-supernovae stage the evolution of the 
specific angular momentum. They find in particular that 
the models with magnetic field lead to neutron stars spin- 
ning about an order of magnitude slower than the models 
without magnetic field by Heger et al. H2000I The model 
with magnetic field is still rotating faster than observed in 
young pulsars, while it is to slow for making the coUapsar 
model possible fWooslev MacFadyen et al. IMIT|l . 

The effect of magnetic field in the core collapse have been 
studied by Fryer & Warren 1)2004(1 . who find that mag- 
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Fig. 4. Evolution of the rotation velocities at the surface 
of the various models with an initial mass of 15 M0 with 
Z = 0.02, during the MS phase. The heavy broken line is 
for model V3 the light broken-line for model Ml and the 
continuous line for model M5. 

netic fields do not play a dominant role in the supernova 
explosion mechanism. 

In paper I of this series (Maeder & Mevnet f2003|l . we 
have calculated a model of a 15 M0 with solar composi- 
tion and an initial rotation velocity Vini of 300 km s^^. 
At a stage near the middle of the MS phase, we have ex- 
amined what happens if we "turn on" the magnetic field 
and apply Spruit's expressions (Spruit 2002). The general 
result was that both the transports of angular momentum 
and of chemical elements by Tayler magnetic instability 
was much larger than the transports by meridional circu- 
lation and shear instability. Thus, this leads us to suspect 
that magnetic field is an important component of stellar 
evolution. Here, we want to compare Spruit's solutions 
and the more general developments expressions developed 
here. Also, we want to study the effects of the field through 
the evolution. For that purpose, we follow the MS evolu- 
tion of 4 different types of models of 15 Mq with solar 
composition. These models are labeled as follows: 

— VO: a model with no rotation and no magnetic field. 

— V3: a model with Viai = 300 km-s~^, and no mag- 
netic field, with the physical assumptions adopted by 
Meynet & Maeder [2000) except for the overshooting 
parameter which was taken here equal to 0. 

— Ml: a model with v\ai of 300 km-s~^ and with magnetic 
field calculated according to the expressions given by 
Spruit H2002|l . i.e., for Case when N^^ > Nt and for 
Case 1 when Nfj, < Nt. In Case 1, when the account 
of radiative losses does not increase the magnetic dif- 



fusivity, another case called IP as given by Eq. H31|l is 
considered (Spruit ,20021 ). 

— M5: a model with of 300 km-s~^ and with magnetic 
field according to the expressions of the present paper. 

In both models Ml and M5, we did not include the effects 
of meridional circulation and horizontal turbulent diffu- 
sion, since these effects were found small with respect to 
those due to the magnetic field, (see remarks in Sect. 6). 
In these first basic tests, overshooting and semiconvection 
were not included. However, we keep the effects of shear 
mixing and the hydrostatic effects of rotation, which dis- 
tort the equipotentials and modify the stellar shape, which 
has great consequences on the T^ff of rotating models. To 
make the comparisons more relevant, we do not include in 
model M5 the energy condition discussed in paper I, since 
this condition was not included by Spruit H2002|l . Some 
further theoretical works are still needed about this con- 
dition and this point will be examined in a further work. 

5.2. Evolution of the internal and superficial rotation 

We examine the evolution of the internal profile of fl{r) 
in the various models. Figs. [5] and El show these pro- 
files at various evolutionary stages during the evolution of 
models V3, Ml and M5 respectively. In model V3, we no- 
tice the smooth, but significant decline of Q{r) outside the 
convective core, where rotation is homogeneous. Such an 
evolution of rotation has been studied in detail by Meynet 
& Maeder PTUHjl and Maeder & Meynet lf2D0I|) at solar 
composition and at lower metallicity. The angular velocity 
is first decreasing in the core as mass loss removes angular 
momentum at stellar surface, because enough coupling of 
rotation is ensured by shear mixing and meridional circu- 
lation to largely compensate the acceleration of rotation 
which would result from the moderate core contraction. 
Only at the end of the MS phase, central rotation increases 
due to the dominant effect of central contraction. In the 
outer layers, a smooth but significant decrease of rotation 
between the core and the surface is created by the above- 
mentioned processes. 

In Fig.|21and|3| we notice several interesting new facts 
for models with a magnetic field: 

— The models have almost an internal solid body rota- 
tion, with a difference of only a few percent between 
the core and the bulk of the envelope, with a short 
transition at their interface. The differential rotation 
parameter q is very small, but still different from zero. 
The relative constancy of is the result of the high 
value of the transport coefficient v by the magnetic 
field (see below Fig. (HJ. We notice that these results 
are in agreement with those of Heger et al. l|2004|l . 
who also find essentially solid body rotation until the 
end of the MS phase. Then in later phases, differential 
rotation appears together with a slowing down of the 
core. 

— There is a general slowing down of rotation during evo- 
lution due to two effects: the mass loss which removes 
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Fig. 5. Internal distribution of the hydrogen mass frac- 
tion X as a function of the Lagrangian mass in the model 
V3 (with rotation only) at various stages of the model 
evolution from top to bottom during the MS-phase. 

angular momentum and the expansion of the external 
layers. The strong internal coupling of rotation trans- 
ports this slowing down to the core. 

— Models Ml and M5, although not strictly identical (see 
Fig.^, exhibit similar Sl-distributions. 

— We notice a small decrease of J7 near the surface. This 
is likely resulting from the expansion of the superficial 
layers after the removal of some mass by stellar winds 
at the surface. 

Fig. 01 shows the evolution of the surface velocities for 
models V3, Ml and M5. Wc note the small differences 
between models Ml and M5, showing that the general 
treatment is not identical to Spruit's solutions. We see 
that the fast initial decline of model V3, due to the rapid 
adjustment of an equilibrium profile by meridional circu- 
lation (cf. Meynet & Maeder I2(J(J(J|I is absent in models 
with magnetic field which rotate nearly like a solid body 
from the ZAMS onwards. Thus models Ml and M5 keep 
an almost constant surface velocity, while model V3 also 
has a slight decline during the rest of the MS phase. Very 
likely this behaviour is different for different initial masses, 
because of the differences in the mass loss rates. Thus, a 
complete grid has to be made in order to establish a basis 
for comparisons with the observations. 

5.3. Internal distribution of hydrogen 

The evolution of the internal distribution of hydrogen 
mass fraction is illustrated in Fig. for model V3 and 
in Fig. ini for model M5 with magnetic field, (models Ml 




Fig. 6. Internal distribution of the hydrogen mass fraction 
AT as a function of the Lagrangian mass in the model M5 
(with magnetic field calculated according to expressions of 
the present paper) at various stages of the model evolution 
from top to bottom during the MS-phase. 

and M5 gives very similar curves with only minor differ- 
ences as illustrated by the evolution of the core sizes in 
Fig.EI). We see that at the end of the MS phase model M5 
produces a convective core which is slightly smaller than 
for model V3. The H-distribution outside the core is much 
smoother in model V3: there is more helium outside the 
formal core of this model than in model M5. The reason is 
that in the magnetic models the transport of chemical ele- 
ments is strongly inhibited by the /x-gradient just outside 
the core: the diffusion coefficient for chemical elements de- 
pends on iVp with a power —6, as shown by Eq. H21|l . In 
the models with only rotation, the inhibition by for the 
shear diffusion goes like a power —2. The fact that more 
helium is burnt in models V3 has some consequences for 
the shape of the tracks in the HR diagram and also for 
the lifetimes as shown in Sect. 5.6 below. 

5.4. The Brunt-Vaisala and Alfven frequencies 

Fig. shows the terms and iV^ of the Brunt-Vaisala 
frequency for the model M5 near the middle of the MS 
phase, (the general behaviour is similar for model Ml, 
thus we do not show it). We see that the thermal com- 
ponent dominates through most of the envelope while 
the mean molecular weight component N'^ is much more 
important close to the core. This large iV^ inhibits the 
chemical mixing. 

The ratio lua/Q, lies everywhere between 10^'^ and 
10~^. The magnetic field B^p is nearly constant through 
the star. Its amplitude is about logB^ = 4.4 through- 
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Fig. 7. Internal distribution of the partial Brunt-Vaisala 
frequencies and iV^ for the model M5 at an age of 
8.384 • 10^ yr with central H-content given by Xc — 0.30. 



out with average deviations smaller than about 0.1 dex, 
where the field is expressed in Gauss. As already sug- 
gested by Maeder & Meynet l|2003|l . the field currently due 
to Tayler-Spruit dynamo is rather large in massive stars. 
The situation is different here from the above work. The 
new point here is that even the equilibrium field is large, 
i.e., the field created by the differential rotation which is 
itself strongly reduced by the field in a consistent feed- 
back. The magnetic field decreases fastly near the stellar 
surface, because differential rotation becomes negligible. 
What is left of the field in the photosphere is difficult 
to ascertain without more detailed models. Also, we re- 
call that for consistency with Spruit's equations (|2)02J, 
we have not accounted here for the energy condition pro- 
posed in paper I (Maeder & Mevnet l2003|l . which tends 
to reduce or suppress the field in the outer layers. 

5.5. DifFusion coefficients for the transports of the 
angular momentum and chemical elements 

Fig. IHI shows the coefficients of transport in model M5. 
The continuous, higher curves represent the coefficient 
v for the transport of the angular momentum as given 
by Eq. H45|l . In the outer layers, the value of ly reaches 
about 10^^ cm^- s""'^, which differs by less than one or- 
der of magnitude from the typical values of the diffusion 
coefficient in a convective zone. In the region where the 
/i-gradient becomes significant, the diffusion is reduced to 
1/ « 10^ cm^-s^^, which is still large. These values of 
implies a strong coupling of the rotation motions, which 
is responsible for the near constancy of the angular veloc- 



Fig. 8. Internal distribution of the coefficients of diffusion. 
The coefficient v (Eg. 145(1 for the transport of angular mo- 
mentum by the magnetic field is shown by the continuous 
line. The coefficient of transport of the chemical elements, 
i.e., the magnetic diffusivity 77 given by Eq. (|19() is shown 
by a short-dashed line. The coefficient Dghcar of the shear 
turbulent mixing is shown in comparison as a dotted line. 
The long-dashed line shows the thermal conductivity K. 
The ages of the models are the same as in the model of 
Fig.lZI 

ity fl through the whole stellar interior. As illustrated by 
Eq. gSJ and (gHl, this coefficient behaves like (il/A^^)** in 
this region, which explains the observed decline towards 
the interior. This decline is also favoured by the fact that 

also decreases towards the deep interior as illustrated 
by Fig. [7| We see that v and the thermal diffusivity K are 
of the same order of magnitude in the region above the 
core, while v is larger in the outer layers. This behaviour 
is consistent with the comparison (cf. Table 1 below) of 
the velocities of meridional circulation and of the growth 
of the magnetic instability. 

The coefficient of magnetic diffusivity 77, which acts 
both for the diffusion of the magnetic field as well as for the 
transport of the chemical elements as discussed in Sect. 3, 
is shown by the dashed-line in Fig. |S1 Close to the core, 
the dependence of rjo (Eq. 121(1 goes with the power of —6 
on Nfj^. This explains why in the deep interior, near the 
core, this coefficient becomes very low, being less than 
10'^ cm^- s~^. These regions where r] is very small tend to 
inhibit the transfer of the chemical elements towards the 
surface The coefficient 77 grows rapidly in the very outer 
layers, rising up to 10^^ cm^- s~^. We may check here the 
very small values of the ratio mentioned above. 

A curve with coefficient -Dshoar due to shear turbulence 
is also shown in Fig. |S1 As we have seen, the fi-gradient 
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Fig. 9. The evolution of the mass of the convective core. 
The dotted-hne apphes to model VO, the long-broken (up- 
per) line to model V3, the short-broken line to model Ml, 
the continuous line to model M5. 

is not strictly zero. The value of I?siicar is negligible in 
the outer layers, while it is larger than the coefficient for 
magnetic transport rj in the layers close to the core. The 
reason is that the dependence in N^^ is weaker for -Dshcar 
than for 77, (being in N^'^ for the first and in iV^^ for 
the second). This means that both shear mixing and mag- 
netic effects have to be taken into account. It is clear that 
the interaction of shear instability and magnetic field may 
be rather complex physically with some important feed- 
back effects, particularly in regions where -Dshcar and rj are 
of the same order. The present work does not solve this 
point, but Fig. |H| shows the orders of magnitude of the 
main parameters. The main conclusion is that magnetic 
coupling dominates the dynamics of rotation. This con- 
clusion is in agreement with that by Heger et al. f 2004') , 
who also find an essentially solid body during MS evolu- 
tion. We also notice that from a different physical dynamo 
model MacDonald & MuUan (J2004| suggest that shear in- 
stability is present over most of the interior, a result which 
possibly supports the simultaneous consideration of shear 
transport and magnetic field. 

5.6. Evolution of the core, lifetimes and tracks in the 
HR diagram 

Fig. El shows the evolution of the mass of the convective 
core. As well known, rotation produces an increase of the 
core mass fraction and of the MS lifetime. Here, the MS 
lifetime of model V3 is about 20 % larger than for model 
VO. The models Ml and M5 with both rotation and mag- 
netic field have a core size and a MS lifetime which are in- 
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Fig. 10. Evolutionary tracks in the HR diagram. The 
dotted line refers to model VO, the long-broken line to 
model V3, the short-broken line to model Ml and the 
continuous line to model M5. 

termediate between those of models VO and V3, but closer 
to the non-rotating model VO. This is consistent with the 
results of Figs. [S] and (HJ which show that there is more 
H burnt in models with rotation only. The mixing of the 
chemical elements by Tayler instability increases only very 
slightly the fuel reservoir with respect to a model without 
rotation, since the strong inhibition (Eq. I19|) by the /z- 
gradient limits very much the importance of the mixing. 
Again, we notice some small differences between models 
Ml and M5. The core is larger in Ml than in M5. The 
reason is likely due to the slightly larger values of Ushcar 
in model Ml close to the core. The deviation of model M5 
from the non-rotating case is very small: the lifetime is 
increased by only 2.4 %. 

Fig. ^1 shows the evolutionary tracks in the HR di- 
agram. In interpreting these tracks, one must remember 
that the effects of the distortion of the stellar surface by 
rotation have been accounted for in the definition of the 
Teff. (In models with rotation, we define the average ef- 
fective temperature as the stellar luminosity divided by 
the real distorted surface of the rotating star, cf. Meynet 
& Maeder Il997|l . This definition is equivalent to seeing 
the star with an average orientation, i.e., with an angle 
i between the rotation axis and the line of sight equal to 
the root of the second Legendre polynomial, i.e., about 54 
degrees. This is why all models with rotation start their 
evolution at a place in the HR diagram, which is different 
from that of the model VO with no rotation. 

Model V3 experiences a larger increase of its luminos- 
ity than model VO, as a result of the more extended mixing 
and larger core as shown in Fig.El The growth of luminos- 
ity for models Ml and M5 is smaller than for model V3, 
because mixing is less extended (cf. Fig. ISj). Models Ml 
and M5 experience about the same increase in luminosity 
as model VO. The end of the MS in these two models is 
a bit more shifted to the red by the atmospheric distor- 
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Fig. 11. Time evolution of the surface helium content Yg 
in mass fraction, of the N/0, N/H and N/C in mass frac- 
tion for the various models VO (dotted hne), V3 (long- 
broken line). Ml (short-broken line) and M5 (continuous 
hne) . 



tion than the zero-age models, since the ratio ^ increases 
from about 0.66 to 0.82 during MS evolution. 

We see a difference of the tracks for models Ml and 
M5. A rough agreement is achieved with the limiting cases 
by Spruit H2U02fl . but the differences are not entierly neg- 
ligible. We also notice the "stair-case" track of the mag- 
netic models, particularly for model Ml. This is the result 
of the rather particular evolution of the core mass of this 
model as illustrated in Fig. |H1 initially model Ml evolves 
close to the model V3 with rotation, then it behaves more 
as the model VO with no rotation, account being given to 
the atmospheric effect. This behaviour is rather different 
for model M5 which evolves like the model without rota- 
tion. This illustrates the different results obtained from 
the asymptotic relations and the general solution. 

5. 7. Evolution of the abundances of helium and CNO 
elements at the stellar surface 

Fig- E] reflects the evolution of the abundances of helium 
and CNO elements at the surface of the various models. 
For model VO, as is well known the enrichment only occurs 
when convective dredge-up appears in the red-supergiant 
stage. For model V3, the growth of the nitrogen at the stel- 
lar surface amounts to about a factor of 2 during the MS 
phase, as a result of shear mixing mainly and of meridional 
circulation to a smaller extent. The growth of the helium 
content is quite small and amounts to about 0.004 in Y^. 
Model Ml shows a late N-enrichment up to the level of 
model V3 near the end of the MS, while it is negligible 
before. The helium enhancement is even higher than for 
model V3. The striking point is that model M5 shows no 
enrichment at all. This is likely a consequence of slightly 
lower Dshear close to the core than for model Ml. There 



the difference between models Ml and M5 is rather critical 
in view of the observational tests to be performed. 

We know from spectroscopic observations that signifi- 
cant surface enrichments in helium and nitrogen occur in 
rotating 0-stars as well in some early B-type stars (cf. 
Howarth & SmithlDUT Villamariz et al.'2002). Thus, the 
fact that the rotating model M5 does not produce such 
enrichments is a negative indication, suggesting that the 
present models with magnetic field do not well correspond 
to reality. The fact that the agreement is better realized 
with model Ml is not meaningful, since this results from 
the use of asymptotic relations which lead to an overesti- 
mate of the diffusion coefficients in some parts of the star. 
Incidentally, this also shows that wrong conclusions could 
be drawn, if the general expressions are not used. The 
models with only rotation much better agree with obser- 
vations (cf. for example Maeder & Mevnet F2001|l . However, 
one must be careful. These comparisons do not necessar- 
ily exclude the presence of any magnetic fields playing a 
role in massive stars. In that respect, the remarks made in 
Sect. 6 below on the equilibrium reached between merid- 
ional circulation and magnetic transport certainly leads to 
more differential rotation and mixing than in the present 
models. 

6. Discussion: interaction between Spruit's 
dynamo and meridional circulation. 

The aim of this section is to discuss the role of the mag- 
netic field with respect to other instabilities, in particu- 
lar the thermal instability, a point which has generally 
not been considered. The various authors, for example 
Heger et al. I|2004|l . usually just sum up the various ef- 
fects. This simple approach may nevertheless have the 
interest to show the order of magnitude of the effects 
considered. When differences of several orders of mag- 
nitude appear between, for example, the coefficients for 
the transport of angular momentum v and -Dshcar, we 
may likely reach a conclusion about which effect domi- 
nates. In this context, we recall that Charbonneau l)2001|l 
has examined the delicate interplay between dynamo and 
meridional circulation. With a different dynamo model, 
he has concluded that the dynamo action remains proba- 
ble in presence of meridional circulation. Our conclusion 
below is not different. If only a dynamo is present, the 
star reaches solid body rotation, which would drive a cir- 
culation. On the contrary, we will see that, if we ignore 
the dynamo, the meridional circulation and star evolution 
build a differential rotation sufficient to make the dynamo 
active. Therefore, we may conclude that in the time evo- 
lution both effects of dynamo and meridional circulation 
are present, with a delicate balance between them. This 
first approach nevertheless still leaves room for a thor- 
ough analysis of the magnetic and thermal instabilities at 
a fundamental level. 

We have seen in paper I (Maeder & Meynet I2003|l 
that, if in the course of evolutionary models calculated 
with only the current effects of rotation (hydrostatic ef- 
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fects, meridional circulation, shear mixing, horizontal tur- 
bulence, etc..) and which have a significant amount of dif- 
ferential rotation, we estimate the magnetic field created 
by Tayler-Spruit dynamo, then the magnetic field and its 
transport effects are found to be very large. The transport 
coefficients v and r\ by the magnetic field were found to 
overcome the above mentioned effects of rotation, in par- 
ticular meridional circulation which is usually the most 
important process for the transport of angular momen- 
tum. This lead us to suggest that magnetic fields may be 
dominating. The situation is probably not so simple. 

In Sect. 5 above, we have followed the evolution of a 
15 Mq star with the magnetic field created by the Tayler- 
Spruit dynamo and we have found an internal equilibrium 
of rotation characterized by a nearly constant 11 (r) in the 
stellar interior due to the strong magnetic coupling ex- 
pressed by the coefficient v. We may now wonder whether 
in these models with « const, the effects of the mag- 
netic field are always larger than those which could be due 
to the thermal instability driving meridional circulation. 
If the answer is "yes" , then clearly Tayler-instability and 
magnetic fields largely dominate the thermal instability 
responsible for meridional circulation. If on the contrary 
the answer is "no" , this means that the situation is more 
complicated. 

The above worry is quite justified, since we know that 
there are various feedback mechanisms playing a role in 
the problem. Models with only rotation lead to a strong 
differential rotation, which in turn may feed the dynamo 
mechanism and the build-up of a magnetic field, which 
then would tend to dampen the differential rotation, 
which is its source. We also know from previous work 
(Meynet & Maeder I2()()()jl that if rotation is constant in 
the interior, as is the case in the initial stellar models 
on the ZAMS, then the thermal instability and the 
resulting velocity U{r) of meridional circulation is much 
larger than in the case for the equilibrium profile i~l{r) 
of rotation. Moreover U (r) is positive everywhere in 
the envelope. This means that angular momentum is 
transported from the outer regions of the star to the 
inner ones, counterbalancing the effects of the transport 
by diffusion. Thus this large meridional circulation will 
enhance differential rotation. The complete feedback loop 
is the following, with two parts A and B: 

A) Build-up of magnetic field and constant rotation: 

differential rotation ^ dynamo => 
magnetic field => il ~ const. 

B) Build-up of circulation and differential rotation: 

Q ^ const. stronger thermal instability ^ 

higher U{r) differential rotation ... 

The basic equation governing the evolution of ^{r) is the 
equation describing the conservation of the angular mo- 



mentum, which is in Lagrangian coordinates (cf. Zahn, 
[BMl Maeder & Zahn lBMIl 



dr 



dr 



(53) 



f2(r) is the average angular velocity on the equipotential. 
If an equilibrium of the above feedback loop is reached in 
a time relatively short with respect to the evolutionary 
timescale, this equation is equal to zero. We would thus 
have a balance of the flux due to meridional circulation 
U{r) by the flux due to the magnetic diffusion expressed 
by V, as shown by Eq. H53(I .The equilibrium is achieved for 



U{r) - 5 - 

r V 9 In r 



(54) 



The factor 5 results from the integration of the angu- 
lar momentum over the isobar, which leads to the above 
Eq. . On the right side of Eq. ifSUl , we have a velocity 
which can be associated with the magnetic diffusion. Let 
us call it t/niagn ^-iid the left member to make the 

distinction. In Table 1, we compare the values of i7magn 
and [/circ at the middle of the MS phase in the model of 
15 M0 studied in Sect. 5. 



Table 1. Comparison between Uc\ 
Mq model at Xc = 0.30. 



\g\ 



and (7magn in a 15 



log r 



8.02 1.97O-02 2.03C+11 0.45c-3 11.1060 3.60c-03 
10.05 1.55C-02 2.24c-(-12 0.62g-4 11.1897 4.46g-03 
13.04 1.85C-02 7.780-I-13 0.49o-4 11.3298 9.00o-02 



We see that Ucu-c dominates over J7,nagn in the deep en- 
velope, while the opposite situation occurs closer to the 
surface; the difference reaches about a factor of 5. Models 
with only rotation as in paper I indicate that the magnetic 
fields should be large, while models with magnetic fields 
as here indicate that the meridional circulation is signifi- 
cant. The basic reason for the difference is the ft{r) profile, 
with important differential rotation in the first case and 
constant rotation in the second. 

We think that an equilibrium profile of il(r) is reached 
as a result of the interaction of the two effects, neither with 
as much differential rotation as in the models with rota- 
tion only (cf. Fig. 1), nor with a constant angular velocity 
as in the present models with magnetic field (cf. Figs. 2 
and 3). This intermediate fl{r) profile will lead to more 
mixing than in the present case. Will this be enough to 
be in agreement with the observations? We do not know 
yet. We emphasize that the equilibrium profile of n{r) is 
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likely defined by a particular balance between J/magn and 
fcirc, with account taken also of the stellar contraction 
and angular momentum conservation. There is probably 
an additional constraint imposed by the fact that the en- 
ergy of the magnetic field originates from the excess energy 
in differential rotation and thus the energy density in the 
magnetic field is limited by its energy source. Such effects 
will be examined in further work. 

7. Conclusions 

The equations for the general case of the Tayler-Spruit dy- 
namo with /z-gradients and nonadiabatic effects have been 
developed. The general expressions consistently allow us 
to recover Spruit's asymptotic relations H2()()2|l . when Nf_i 
or Nt are zero. Numerically the general and asymptotic 
results are consistent with some quantitative differences. 
This demonstrates the need to use the general solutions, 
in particular if we want to perform observational tests. 

Some initial results on the role of possible magnetic 
field in stellar models have been found. The Tayler-Spruit 
dynamo imposes that the stars to rotate nearly as a solid 
body. This leads to a different evolution of rotational ve- 
locities as stars move away from the zero-age sequence. In 
general, the rotation velocities are higher when a magnetic 
field is present. A second general result is that internal 
mixing and surface enrichments in products of the CNO 
burning are smaller when a magnetic field is present than 
in rotating models without a magnetic field. Since spec- 
troscopic observations of OB stars show significant N- and 
He-enrichments, this cast some doubt as to whether mag- 
netic fields have a dominant effect, especially more than 
the models which have only rotation and no magnetic field 
show a good agreement with observations. 

If magnetic fields are present, we may expect, as shown 
in Sect. 5, that the magnetic instability and the thermal 
instability which drives meridional circulation reach some 
balance in their complex feedback loop. For reasons ex- 
plained in Sect. 5, this balance should be characterized 
by less differential rotation than in current models with 
only rotation and by more differential rotation than in 
the present models with a magnetic field. Such models 
may lead to more mixing than the present ones. 

Finally, we emphasize that at this stage we really do 
not know how the Tayler-Spruit instability interacts with 
the meridional circulation and how the energy constraints 
may modify the results. But we have explored some conse- 
quences of the magnetic field in the absence of meridional 
circulation and examined the possibilities for its existence. 
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